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The multipocket Fermi surfaces of iron-based superconductors promote pairing states with both 
s±-wave and (ia,2_^2-wave symmetry. We argue that the competition between these two order 
parameters could lead to a time-reversal-symmetry breaking state with s + id-pairing symmetry in 
the iron-based superconductors, and propose serveral scenarios in which this phase may be found. 
To understand the emergence of such a pairing state on a more rigorous footing, we start from a 
microscopic 5-orbital description representative for the pnictides. Using a combined approach of 
functional renormalization group and mean-field analysis, we identify the microscopic parameters of 
the 8 + i(i-pairing state. There, we find the most promising region for s + icZ-pairing in the electron 
doped regime with an enhanced pnictogen height. 

PACS numbers: 74.20.Mn, 74.20.Rp, 74.25.Jb, 74.72.Jb 



Iron based superconductors (SC) offer an appealing 
platform to investigate the interplay among pairing in- 
teractions, pairing symmetries and Fermi surface topolo- 
gies pj. Generally, repulsive interactions in momentum 
space can lead to a change of sign in the pairing ampli- 
tude. A large class of iron based SC have disconnected 
Fermi surface pockets, consisting of hole pockets at the 
r = (0, 0) and possibly M = (tt, tt) points, and two elec- 
tron pockets at the X = (tt, 0)/(0, tt) points in the un- 
folded Brillouin zone (BZ) with one iron atom per unit 
cell. When the repulsive interactions between the hole 
and the electron pockets dominate, an s± pairing symme- 
try can be obtained [IHJ. On the other hand, when the 
repulsive interactions between the two electron pockets 
dominate, a propensity toward d-wave pairing symmetry 
can be expected. When both types of interactions are 
comparable, there is hence a frustration between the two 
pairing symmetry types. A recent theoretical proposal 
suggests that the system can resolve the frustration by 
a pairing state with the s -\- id pairing symmetry which 
spontaneously breaks time-reversal (TR) symmetry [5]. 
The possibility of a TR-symmetry breaking pairing state 
due to frustrating pairing interactions among three or 
more Fermi pockets has also been investigated in several 
other contexts P^TO]. In general, time reversal breaking 
pairing states have rather accessible experimental signa- 
tures, and several proposals have been suggested in the 
context of iron based SC [5 . 

In principle there are various experimentally tunable 
parameters to drive the competition between s-wave and 
(i-wave in the pnictides, giving the opportunity to start 
from both limits. In Ka^Bai_a;Fe2As2, the Fermi surface 
topology can be chosen as a paradigmatic setup for s±, 
consisting of hole pockets at F and the electron pock- 
ets at X for optimal doping x 0.4. For increasing x. 
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FIG. 1. (Color online). Frustrating the d-wave limits of 
KFe2As2 (a) and KxFe2-ySe2 (b). Upon doping or differ- 
ently induced band structure effects, electron pockets appear 
(dashed red) in (a) and a hole pocket appears (dashed red) 
in (b) which populate the q ^ (tt, 0)/(0, tt) scattering chan- 
nels enhancing the s± symmetry. This leads to frustration 
providing the background for s -\- id pairing. 



however, the electron pockets decrease, and have nearly 
disappeared for x = 1 [Fig. [l], which has been recently 
suggested to host a (i-wave pairing symmetry pjj. In 
this system, it is hence plausible that 8i s -\- id pairing 
state can be realized for intermediate values of x. In the 
chalcogenide KxFe2-ySe2^ the electron pockets at the X 
points dominate, and, for a situation seemingly inverse 
to KFe2As2, a (i-wave pairing symmetry may likewise be 
expected [121 US]- (It should be noted that the actual 
pairing symmetry in the chalcogenides is far from set- 
tled, as a strong coupling perspective may likewise sug- 
gest s± pairing [14 .) By tuning doping or other possible 
parameters affecting the band structure such as pressure, 
one possibly induces a pocket at F, increasing the ten- 
dency towards s± pairing symmetry [Fig.fTl. In this case. 
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FIG. 2. (Color online). Competing pairing orders and SC form factors for Ui{dx'2_Y'^) — — 2.5ey (a) and Ui{dx2_Y'^) — 
1.6U* (b) at electron doped filling n = 6.13. RG channel flow (al,bl) and s±-gap form factor (a2,b2). s±/(i^2_^2 transition 
from (a) to (b): increasing Ui{dx2_Y'^) enhances the gap anisotropy of the s±-form factor on the electron pockets [fc-patching: 
points 33-64 see (a3)] shown in (a2,b2) until its symmetry switches to d^2_y2. (a3,b3) Interactions mediated by ?7i, setting up 
the s±-pairing tendency between hole pockets and electron pockets (F o X) and the competing (i3,2_ ^2 -pairing symmetry due 
to electron-electron {X ^ X) interaction, (c) The variation of the pnictogen height hp particularly affects the spread of the 
(ix2-Y2 -orbital and therefore Ui{dx'2-Y'^)i as it is oriented to the planar projection of the pnictogen. 



one could also expect an s + pairing state. By sys- 
tematically tuning the Fermi pocket topologies, one can 
compare the predicted pairing symmetries with exper- 
iments, and determine the nature of the pairing inter- 
action by these investigations, starting from compound 
settings with a suspected d-wave symmetry [Fig. [l]. 

In the following, we rather intend to start from an 
s± pairing state to begin with, and address how we can 
enhance the competitiveness of the d-wave symmetry to 
drive the system into the s ^ id regime. The reason for 
this is two- fold. First, the s± symmetry is much more 
generic for the different classes of pnictides. Second, as 
we will see below, we find the most promising setup to 
be located on the electron doped side of pnictides, where 



high-quality samples have already been grown for dif- 
ferent families. We hence believe that this regime may 
be the experimentally most accessible scenario at the 
present stage, which is why we explicate it in detail. In 
this paper, we investigate the microscopic mechanism of 
the s-\-id pairing state by the functional renormalization 
group (£RG) method of a five band model. We system- 
atically vary the doping level and the strength of intra- 
orbital interaction, which determine the ratio between 
the electron-hole pocket and the electron-electron pocket 
mediated pairing interactions. In this microscopic inves- 
tigation, we find that the s -\- id pairing state can be 
realized in the intermediate electron-doped regime, given 
that we also adjust the pnictogen height parameter of the 
system appropriately. 



We start from a representative 5-band model for the 
pnictides which is obtained from LDA-type calcula- 
tions [3] . It has been considered by us before as a starting 
point for explaining the difference between the isovalent 
P-based and As-based pnictides [15 . The LDA "non- 



interacting" part is given by 

5 

^0 = E E 4a«^a6(fc)c,,,. (1) 

k,s 0,6=1 

Here c's stand for electron annihilation operators, a, b for 
the (i-orbitals, s denotes the spin indices and Kab{k) the 
orbital (i.e. maximally-localized Wannier function) ma- 
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trix elements of the Kohn-Sham Hamiltonian. The band 
structure features electron pockets at X and hole pock- 
ets at r, which is the typical situation in the pnictides 
[Fig. [2] for sufficient electron doping. The many-body 
interaction part is given by the intra- and inter-orbital 
interactions Ui and /72, as well as the Hund's coupling 
Jh and the pair hopping Jpair- 
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where rii^as denote density operators at site i of spin s 
in orbital a. Typical interaction settings are dominated 
by intra-orbital coupling, Ui > U2 > Ju ^ ^pair- In 
the fRG [T5H2O], one starts from the bare many-body 
interaction Q in the Hamiltonian. The pairing is dy- 
namically generated by systematically integrating out the 
high-energy degrees of freedom including important fluc- 
tuations (magnetic, SC, screening, vertex corrections) on 
equal footing. This differs from the RPA which takes 
right from the outset a magnetically driven SF-type of 
pairing interaction. For a given instability character- 
ized by some order parameter Ok, the effective inter- 
action vertex VA(ki, k2, ka, in the particular order- 
ing channel can be written in shorthand notation as 
p VA(k, p)[O^Op]. Accordingly, the effective inter- 
action vertex VA(k, — k, p, — p)) in the Cooper channel 
can be decomposed into different eigenmode contribu- 
tions [isi nil 

Ff^(k,p) = ^cr(A)/^^'^(k)*/^^'^(p), (3) 

i 

where i is a symmetry decomposition index, and the 
leading instability of that channel corresponds to an 
eigenvalue cf^(A) ffist diverging under the flow of A. 
/^^'*(k) is the SC form factor of pairing mode i which 
tells us about the SC pairing symmetry and hence gap 
structure associated with it. In the fRG, from the fi- 
nal Cooper channel in the effective interaction vertex, 
this quantity is computed along the discretized Fermi 
surfaces [Fig. [2|a3)], and the leading SC instabilities 
are plotted in Figs. |2|al) and (bl). The interaction 
parameters are kept fixed at the representative setup 
Ui = 2.5eV,U2 = 2.2eF, Jh = 1.2eF, Jpair = 0.2eV {Ui 
for the (ix2-y 2 -orbital is varied as explicitly stated in 
Figs. [2] and [3|. The relatively large bare value of Jh 
is motivated partly by recent findings, in particular, for 
a sizable Hund's rule coupling ^fl Furthermore, 
as a parameter trend, larger Jh and smaller Jpair tends 
to prefer the s + zd-phase in the electron-doped regime 
for rather moderate values of intra-orbital coupling Ui 
[Fig. I- 
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FIG. 3. (Color online). Preferred pairing as a function 
of electron doping and intra-orbital Coulomb interaction 
Ui{dx'2_Y'^)- The results are obtained by minimizing the 
mean-field free energy of the effective theory taken from fRG 
at A ^ .OOley. At 27% electron doping, the s + id-pairing 
state occurs at Ui{dx'2_Y2) = SeV which is comparable to the 
intra-orbital repulsion in the remaining orbit als Ui — 2.5eV. 



The situation in Fig. |2] is representative for moderate 
electron doping and interaction scales of the pnictides, 
where the T ^ X pair scattering between the hole pock- 
ets at F and the electron pockets at X dominates. Al- 
ready from the BCS gap equation, a finite momentum 
transfer can induce pairing only when the wave vector 
of such an interaction connects regions on one FS (such 
as in the cuprate case), or regions on different FSs (such 
as in the pnictide case), which have opposite signs of 
the SC order parameter. This corresponds to putting 
the electron pairs in an anisotropic wave function such 
as sign-reversing 5- wave {s±) in Fig. [2|l, where the wave 
vector (tt, 0) in the unfolded BZ connects hole and elec- 
tron pockets with a sign-changing s± gap [2l[4]. However, 
in the fRG calculation of Fig. (2)3 with increased Ui inter- 
action on the (ix2-y2 orbital, a green arrow for X ^ X 
scattering indicates additional interactions that become 
similarly important as the (tt, 0) channel. This increased 
Ui can be tuned by the pnictogen height as explained be- 
low and frustrates the previous "pure" s± limit (F ^ X). 
The system then strikes a compromise [15 , 23 by enhanc- 
ing the anisotropy of the formfactor (denoted by /^^(k) 
in Fig. |2| on the electron pockets at X. The multi-band 
SC hence adjusts the momentum dependence of the gap, 
i.e. its anisotropy, so as to minimize the effect of the 
Coulomb repulsion which is set up by the competition 
between s± and (ia,2_^2-wave channels. 

We now have all ingredients to tune the pairing sym- 
metry from 5±-wave to (ia.2_^2-wave and, eventually, into 
the TR-symmetry broken s -\- id phase. In most of the 
iron-based SC, the tendency towards 5±-pairing occurs 
slightly more pronounced than the competing dx2_y2- 
pairing and, at first glance, the resulting frustration ap- 
pears to be too small for causing s + id-pairing. There- 
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fore, in order to increase frustration, we somehow have to 
enhance the pair-scattering between the electron pock- 
ets at X which then promotes the sub-leading dx2_y2- 
channel. As shown in a-priori determinations of the in- 
teraction in Eq. expressed in terms of orbital ma- 
trix elements, the pnictogen height hp (measured from 
the Fe-plane [Fig. [2]3]) has a substantial influence on the 
intra-orbital interaction Ui between (ix2-y 2 -Wannier or- 
bitals [24] , which can be either modified by isovalent dop- 
ing or pressure. By increasing hp^ the Wannier functions 
in this orbital are further localized, causing an increase 
of Ui{dx'2-Y^)- In Fig- we have already used this 
fact to demonstrate that, for moderate e-doping (13%), 
large values of this matrix element drive the SC instabil- 
ity from s±— to dx2_y2-wdive symmetry. 

For this setup, we present our predictions for TR- 
symmetry breaking in a schematic phase diagram in 
Fig. |3j where we plot the leading s±, dx2-y2 and finally 
s -\- id SC solutions as a function of Ui{dx'2-Y'^)^ and 
electron doping. There, we have used our fRG result 
as a starting point for a renormalized mean field anal- 
ysis |25|. In this MF+fRG approach, the one-loop fiow 
is stopped at a scale A which is small compared to the 
bandwidth, but still safely above the scale Ac, where the 
2-particle vertex diverges. In this range, the particular 
choice of the cutoff A does not significantly infiuence the 
results in Fig. [3j The renormalized coupling function 
V^{ki, ^2, ^3, ^4) is taken as an input for the mean field 
treatment of the remaining modes. As shown in Fig. |2j 
the regime of 5±/(i-wave pairing competition features a 
single channel SC instability without other competing 
(e.g. magnetic) instabilities and, therefore, justifies 

V'^ik^MMM) - VP^'^{kuks)Sk,,-k,Sk,,-ks. (4) 

with V^^^^(/ci,/c3) = V^{ki,-ki,ks,-ks). The ef- 
fective theory for quasiparticles near the Fermi surface 
(|^(/c)| < A) is modeled by the reduced Hamiltonian 



ks 



k,q 



(5) 

where ^{k) is taken as the bare dispersion due to only 
weak band-renormalization effects. The MF solution of 
this reduced Hamiltonian is obtained as in BCS theory, 
by solving the self-consistent gap-equation and calculat- 
ing the corresponding grand potential which is 



^stat ^ 



|Afcp + 2e(/c)2 



^2V^(/c)2 + |AfcP 



Within a reasonable range of parameters for the electron- 
doped pnictides, we then find a regime favoring s + id- 
pairing due to 



Note, that this phase regime is only a lower bound for 
the existence of s -\- id which probably is much larger. 
This is because the fRG setup at present only allows us 
to obtain the leading SC instability at some finite Ac, 
while the s -\- id phase may well set in below Ac. This 
would manifest itself as a change of the SC phase as a 
function of temperature. 

In summary, we have presented a microscopic analy- 
sis, based on a-priori electronic structure determinations 
and a combination of the fRG with an MF treatment 
of the remaining low-energy states, to derive a kind of 
"guiding principle" for a possible s -\- id pairing state in 
the pnictides. For the case of increased electron doping 
and pnictogen height, we have illustrated how this drives 
the system into an s ^ id SC state. Aside from this ex- 
ample, other regimes in the pnictides likewise promise 
the possible realization of an s + id state, such as hole- 
doped (K,Ba)-122 interpolating between the 5- wave limit 
{x ~ 0.4) and d-wave limit (x ~ 1) as well as possibly 
the chalcogenides KxFe2-ySe2' 
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